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ABSTRACT 
An  approximate  method  is  presented  which  yields  the  supersonic 
nozzle  profile  necessary  to  keep  an  expanding  high  temperature  gas  mix- 
ture approximately  in  chemical  equilibriurrio     From  equilibrium  calculations 
for  given  chamber  conditions,   the  freeze  area  ratio     is  determined  by  use 
of  a  modified  Bray  criterion.    A  plot  of  area  ratio  versus  a  function  of  noz- 
zle geometry  defines  a  limiting  nozzle  angle  beyond  which  freezing  is 
forecast  to  occur.     This  permits fen,  ,expa*s ion  angle  at  each  area  ratio 
to  be  chosen  to  keep  the  flow  near  equilibrium.     The  effectiveness  of  the 
resulting  nozzle  profile  is  evaluated  by  comparing  exact  calculations 
with  equilibrium  calculations.    Contouring  is  shown  to  be  effective  for 
small  nozzles  when  a  correction  factor  is  applied.     The  hydrogen-oxygen 
system  is  used  as  an  example. 
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1.         Introduction. 

During  the  last  decade  a  marked  increase  of  interest  in  rocket  design 
and  performance  prediction  has  occurred.     Payloads  for  space  ^exploration 
missions  often  differ  greatly  from  those  of  a  military  nature,   but  the  vehicles 
involved  have  a  major  common  factor  which  is  the  rocket  engine.    Methods 
of  propulsion  currently  under  investigation  for  the  future  include  ion,  mag- 
netoplasma,  and  nuclear  fission;  however,  the  primary  type  today  is  the 
chemical  rocket  engine. 

The  idea  of  jet  propulsion  is  old.    In  the  chemical  rocket,   potential 
chemical  energy  of  the  fuel  and  oxidizer  is  converted  into  thermal  energy  in 
the  combustion  chamber.    A  convergent-divergent  de  Laval  nozzle  then  accel- 
erates the  gas  stream  from  subsonic  to  supersonic  velocity.    This  converts 
the  thermal  energy  and  pressure  potential  energy  of  the  nearly  stagnant,  hot 
gas  ahead  of  the  nozzle  into  the  kinetic  energy  of  the  supersonic  stream 
which  emerges  with  lower  temperature  and  pressure.    A  reaction  is  imparted 
which  is  directly  related  to  the  time  rate  of  change  of  the  momentum  of  the,. 
ejected  matter. 

At  the  high  temperatures  encountered  in  modem  combustion  chambers, 
a  large  part  of  the  chemical  energy  released  may  be  used  in  breaking  chemical 
bonds  to  produce  simpler  species  (dissociation).    As  the  temperature  and 
pressure  decrease  during  expansion,   recombination  of  the  simpler  species 
occurs  with  a  resulting  release  of  chemically  bound  energy  to  the  gas  flow0 
A  significant  portion  of  the  chemical  energy  is  not  available  for  conversion 
to  kinetic  energy  unless  reassociation  of  the  dissociated  combustion  products 
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occur.  Thus,  the  thrust  produced  by  any  propulsive  system  in  which  combus- 
tion products  are  expanded  through  a  nozzle  will  be  reduced  if  the  recombinar- 
tion  process  cannot  occur  rapidly  enough  to  remain  in  equilibrium. 

In  a  rocket  nozzle  the  elementary  reactions  between  the  dissociated 
combustion  products  seldom  proceed  at  rates  which  are  sufficient  to  main- 
tain true  chemical  equilibrium  throughout  the  nozzle  expansion  £l6j*.     For 
a  hypersonic  ramjet,   calculations  indicate  net  thrusts  which  are  attractively 
high  if  equilibrium  is  maintained,   but  which  approach  zero  around  mach  num- 
bers of  ten  or  eleven  if  the  combustion  products  are  frozen,    i.e. ,   constant 
composition  £3jf . 

To  facilitate  further  discussion,  the  following  descriptions  of  flows  are 
given: 

Frozen:  The  chemical  reaction  rates  are  zero  with  composition 

invariant.     Typical  non-reacting  flows  come  under  this 
classification. 
Equilibrium:   The  chemical  reaction  rates  are  infinite,  and  composition 
of  the  gas  adjusts  instantaneously  to  changes  in  temper- 
ature and  pressure. 
Exact:  The  chemical  reaction  rates  are  finite.    A  finite  time  is 

required  for  the  composition  of  the  gas  to  adjust  to  the 
changes  dictated  by  the  pressure  and  temperature  variations 
in  the  nozzle. 

*  Numbers  in  brackets  indicate  references  on  page30u 
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Near-equilibrium  flow  and  near-frozen  flow  imply  a  small  deviation 
from  the  artificial  cases  of  zero  or  infinite  reaction  rates  without  an  appreci- 
able  change  in  the  resulting  performance  prediction.     Equilibrium  flow  yields 
maximum  performance  while  frozen  flow  yields  minimum  performance. 

Consider  the  following  typical  three-body  reaction: 

A 

A    +    B   +    M      ^=      AB    +    M 
R" 

where  M  is  a  third  body  which  removes  the  energy  during  recombination,  and 

R°  and  R"  are  the  rates  of  the  reaction  in  the  forward  (recombination)  and 

backward  (dissociation)  directions  respectively.     In  the  combustion  chamber, 

the  reaction  would  move  to  the  left.    At  lower  temperatures  and  pressures,  the 

equilibrium  composition  would  shift  toward  the  formation  of  AB„ 

If,  however,  the  reaction  rates  are  not  fast  enough  to  maintain  the 
equilibrium  conditions,   recombination  will  lag.     The  gas  leaving  the  nozzle 
will  have  some  dissociated  products  which  have  not  contributed  to  the  kinetic 
energy  of  the  flow.    A  reduction  in  thrust  will  occur  if  the  recombination  pro- 
cess does  not  occur  rapidly  enough  to  remain  in  equilibrium. 

It  is  assumed  that  the  working  fluid  is  a  mixture  of  ideal  gases;  that 
the  flow  is  steady,  adiabatic,  and  quasi-one-dimensional;  and  that  the  effects 
of  diffusion,   heat  conduction,  and  viscosity  are  negligible.     Equilibrium  flow 
calculations  involve  the  simultaneous  solution  of  the  algebraic  equations 
for  conservation  of  mass,   energy,   momentum,  and  species  plus  the  equation 
of  state.    Frozen  flow  calculations  are  the  same  except  that  the  species 
composition  is  constant  during  the  flow.    These  two  types  of  calculations 
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present  no  particular  difficulties.     Many  computer  programs  have  been  written 
to  uyiieid  gas  dynamic  and  thermodjjmamic  data  at  selected  stations  |~4„  5,  6}  . 

For  an  exact  flow,  though,  the  problem  is  more  complex.    Sarli,  etaL, , 
[16]  ,   summarizes  the  difficulties.    The  calculations  require  the  simultaneous 
solution  of  the  mass,   momentum,  and  energy  conservation  equations;  state 
and  mixture  enthalpy  equations;  kinetic  differential  equations  for  the  finite 
rate  of  production  of  molecular  species;  and  atomic  continuity  equations  for 
all  atoms  present.    The  equations  used  in  the  exact  solution  are  given  in 
Appendix  I.    Again  the  assumptions  are  that  the  flow  is  one-dimensional, 
inviscid,  and  adiabatic,  with. negligible  diffusion.     The  simultaneous  solu- 
tion of  the  above  equations  requires  numerical  techniques;  hence,  a  digital 
computer  is  desirable  [9,  1.0,  16}.  . 

Three  basic  problems  exist,  however,  regarding  these  exact  calculations 
[l6]0    These  are:    1.     The  numerical  integration  of  the  equations,  when  near- 
equilibrium  conditions  exist  in  the  flow,  is  relatively  slow  and  potentially 
inaccurate;    2.      It  is  necessary  to  specify  the  detailed  mechanism  and  in- 
dividual rates  for  the  reactions  taking  place;    3.      The  cost  of  the  computer 
time  involved  prohibits  casual  investigations  of  complex,   reacting  flows. 

Today's  space  ventures  require  thrust  prediction  and/or  weight  minimi- 
zation such  that  a  design  based  solely  on  a  frozen  flow  calculation  or  an 
equilibrium  flow  calculation  would  be  unacceptable.     This  does  not  mean 
that  equilibrium  and  frozen  calculations  are  to  be  abandoned,,     In  view  of  the 
computational  problems  mentioned  above  that  exist  with  the  exact  solutions, 
a  method  of  establishing  guidelines  for  the  design  of  rocket  nozzles  is  desir- 
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able  which  utilizes  approximate  methods. 

One  of  the  first  to  investigate  a  modified  equilibrium/frozen  calcula- 
tion to  predict  actual  performance  was  Penner  [12].    Criteria  were  established 
which,   in  some  instances,   permitted  a  nozzle  flow  to  be  classified  as  near- 
equilibrium  or  near-frozen.     Thus,  for  a  given  set  of  conditions,  a  reasonable 
prediction  of  performance  could  be  given  without  encountering  the  complexion 
tiesoffan  exact  calculation.      Reichenbach   Il5]  proposed  an  iterative  proced- 
ure for  the  solution  of  nozzle  flow  problems.    Agreement  with  exact  calcula- 
tions was  generally  good.     Several  investigators  1 1,    ll]  have  considered  the 
design  of  an  optimum  nozzle. 

Another  early  investigator  in  this  area  was  Bray  [~2,   3J.     His  work,  has 
received  considerable  attention  since  its  first  presentation  in  1958.     Many 
of  the  current  papers  in  the  field  of  nozzle  flow  with  chemical  reactions 
utilize  his  concepts. 

The  Bray  criterion  is  an  approximate  procedure  for  prediciting  where  a 
reaction  has  departed  significantly  from  equilibrium  in  a  reacting  nozzle  flow. 
This  prediction  is  made  by  determining  the  point  at  which  the  forward  rate  of 
reaction,   R",   becomes  of  the  same  order  as  the  rate  required  to  maintain 
equilibrium.     The  flow  upstream  is  considered  in  equilibrium,  while  that  down- 
stream is  considered  frozen  (invariant  in  chemical  composition).    While  Bray's 
original  work,  dealt  with  a  single  reaction  system,  later  investigators  |"3,   7, 
8,    1 6 J  ^extended  the  principle  to  include  multiple  reaction  systems. 

An  important  point  to  consider  is  the  effect  of  nozzle  size  and  geometry 
on  performance.     Penner  points  out  [12]  that  the  reaction  time  must  be  small 
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\compared  with  the  residence  time  to  achieve  near-equilibrium.     Hence,   for 
small  nozzles  (short  residence  time)  it  is  difficult  to  achieve  near-equilibrium 
flows. 

After  the  above  mentioned  calculations  are  completed,   it  is  then  appa- 
rent whether  or  not  the  boundary  configuration  chosen  (nozzle  geometry)  is 
likely  to  give  the  desired  performance,   i.e.  thrust  and  specific  impulse.    A 
series  of  computations  using  different  nozzle  geometries  would  eventually 
yield  a  suitable  nozzle  design,   but  this  would  not  be  economical  in  terms  of 
expense  and  labor. 

Kushida  (8j  and  Koppang,   et  al.,  (7 J  presented  a  method  utilizing  the 
Bray  criterion  to  determine  the  freeze  point  for  a  given  nozzle.    Reichenbach 
[14]  suggested  to  this  investigator  that  a  reverse  approach  might  be  plausible. 
Since  an  estimate  is  available  from  Bray's  criterion  for  the  location  of  a  freeze 
point  as  a  function  of  area  ratio,  a  nozzle  designed  to  avoid  the  conditions 
that  lead  to  freezing  of  the  flow  could  approach  the  optimum  performance  for 
a  given  system.    A  further  advantage  would  be  the  relatively  simple  calcula- 
tions involved  which  use  equilibrium  data. 

An  analysis  of  this  approach  was  undertaken  in  this  paper.     The  system 
investigated  was  hydrogen  and  oxygen.     To  compare  and  evaluate  the  results, 
an  exact  calculation  of  the  system  was  performed. 

This  paper  was  prepared  during  the  1965-1966  academic  year  at  the 
United  States  Naval  Postgraduate  School,  Monterey,  California.   Acknowledge- 
ment is  gratefully  made  to  Dr.  R.  E.  Reichenbach  of  the  Aeronautics  Department 
for?. his  guidance  and  consultation  as  thesis  advisor,  and  to  the  staff  of  the 
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computer  facility  for  their  cooperation  and  assistance  in  adapting  the  com- 
puter programs  to  the  Control  Data  Corporation  1604  system. 
2.        Analysis. 

A  common  form  [V,    8,    16,    17]  for  Bray's  sudden  freezing  criterion  is 
the  following  semi-theoretical  relaxation  parameter: 

Equilibrium  rate  of  change  of  species  1 

Rate  of  the  recombination  reaction  Rj  **' 

Thus,   B  relates  equilibrium  composition  change  to  the  exact  composi- 
tion change.    If  the  ratio  is  large,  the  reaction  would  be  near-frozen;  if  the 
ratio  is  small,  the  reaction  would  be  near-equilibrium.    It  is  generally  agreed 
that  the  freeze  point  is  determined  when  B  is  near  unity  [l,   7,   8j. 

For  a  three-body  recombination  reaction  as  shown  in  equation  (l)„  the 
recombination  rate,  R\  and  the  dissociation  rate,  R",  are  given  by  J7] ; 
R'    =    kfn3XAXBXM  (2) 

R"    =    kbn2xABXM  (3) 

where  kf  and  k^  are  the  reaction  rate  constants  (assumed  to  be  a  function  of 

temperature  only);  n  is  the  molar  density  of  the  gaseous  mixture;  and  X^,  Xg, 
XAB'   anc*  XM  are  ^e  m°le  fractions  of  the  respective  atoms  and  molecules. 
It  has  been  shown  [7 ,   8]  that  the  bimolecular  atom  transfer,  or  meta- 
thetic  reactions,  are  an  order  of  magnitude  faster  than  the  three-body  react- 
ions.   Also,  the  energy  released  by  the  metathetic  reactions  is  small  com- 
pared with  the  three-body  reactions.    Consequently,  the  metathetic  reactions 
are  assumed  to  be  in  -  equilibrium  for  this  study.,    Equation  (4)  is  a  typical 
metathetic  reaction. 
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OH+    0        -.  "     H  +    02  (4) 

R" 

The  main  reactions  to  consider  in  the  hydrogen -oxygen  system  are  the 

following  [7 ,    8,    16]  : 

1.  H+H+M    =    H2+M 

2.  0    +    0    +    M    -    02    +    M  (5) 

3.  H    +   0+M    =    OH+M 

4„       OH    +    H    +    M    =    H20  +  M 
It  has  been  found  [7]  that  a  system  of  this  type  is  dependent  on  the 
dissappearance  of  hydrogen  to  determine  the  freeze  point.    As  shown  by  'h 
Kushida  [8]  and  Koppang,   et  al.  [7],  the  net  rate  of  disappearance  of  monatomic 

hydrogen  is  given  by: 

d(X„/ffi) 

-mn    —f =    R8!  +  Ra3  +  R'4  -  R"  x  -  R"3  -  R"4  (6) 

dt 

where:    m  is  the  mean  molecular  weight,  and  the  subscripts  on  R  refer  to 

the  reactions  in  equation  (5). 

The  reaction  rates  are  determined  as  shown  by  Penner  [l2j. 


R^    =    k^KnXj^CnXj^nX^    =    2  kin3x|x 


M 


R'3    =    k3ll)(nXH)(nX0)(nXM)    =    ^X.^X^  (7) 

R'4    =    k4(l)(nX0^(nX^(nXM)    =    k^X^^ 

Since  Xjyj  is  the  mole  fraction  of  any  third  body,   it  is  unity0     Rewriting 
equation  (7): 

R'      -    2  kj  n3XH2  (8) 
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R>3    "    k3  n    XHX0 
R'4   =    k4  n    XQHXH 

To  obtain  the  expression  for  B  as  shown  in  equation  (l),  the  equilibrium 
rate  of  change  of  hydrogen  concentration  is  divided  by  the  rate  of  the  recom- 
bination reactions  for  hydrogen.     It  should  be  noted  that  R'  represents  the 
rate  of  recombination  and  R"  represents  the  rate  of  dissociation. 


K-i   tK«3+K/ 


=     mn  d  (AH/ni 


2  kln3xH2  +  k3n3xHX0  +  k4n3xOHXH 

=     I  d(xH/m)/dt  (9) 

XHn2(2k1XH+k3X0+k4XOH) 

B   =     dQn  m/m)l  /dt 

n2(2k1XH+k3X0+k4XOH) 

Introducing  v  =  dL/dt  and  dividing  both  sides  by  d(6)  where   €  =  A/A* 

(area  divided  by  the  throat  area): 


d€  n2(2  kxXH  +  k3XQ  +  k4XQH)  d6 

Rearranging: 

B        =    v  dLln  (xH/m)]  /de I 

de/dL         n2  (2  klXH  +  k3XQ  +  k4XQH) 

The  right  side  of  equation  (10)  can  be  calculated  using  equilibrium 

data  for  a  given  system.    The  rate  constants,  kf  and  k^,  are  assumed  to  be 

functions  of  temperature  only.     By  plotting  the  right  side  of  equation  (10) 

versus  area  ratio,  a  means  of  determining  the  freeze  point  is  provided.    For 

every  freeze  area  ratio,  there  is  a  corresponding  value  for  B/(d6'/dL). 

(Curve  A,  Fig.   1) 
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d£/dL 


Equilibrium 
Margin 
(EM)  — 


Freeze  Area  Ratio 
Figure    1:      Freeze  Point  Location 

Since  B  is  assumed  to  be  constant  and  equal  to  unity,   it  is  thus  seen 
that  the  local  value  of  d^/dL  determines  the  freeze  point.    Therefore,  if  the 
nozzle  is  designed  such  that  dt/dL  is  less  than  the  local  value  for  each 
freeze  area  ratio  (Curve  A"),  the  flow  should  remain  in  near-equilibrium.     If 
the  nozzle  is  designed  such  that  d6/dL  is  greater  than  the  local  value  for  each 
freeze  area  ratio  (Curve  A"),  the  flow  should  become  near-frozen. 

The  vertical  distance  between  Curve  A  and  Curve  A'  is  arbitrary  and 
not  necessarily  constant.     It  is  important  to  note,  though,  that  a  small  dis- 
placement from  Curve  A  possibly  would  not  yield  near-equilibrium  flow  be- 
cause of  the  inherent  inaccuracies  introduced  by  the  assumptions.    A  large 
displacement  would  increase  the  nozzle  length  because  of  the  slower  expan- 
sion to  a  given  area  ratio.    This  displacement  will  be  referred  to  as  the 
equilibrium  margin  (EM).     Since  the  curves  in  Fig.    1  are  usually  plotted  on 
semi-logarithmic  paper,  the  constant  displacement  EM  becomes  a  multipler 
to  the  basic  curve. 

The  information  available  from  Fig  1  thus  provides  a  guide  for  nozzle 
design  which  incorporates  the  effects  of  chemistry. 
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3.         Procedure. 

The  first  step  necessary  in  determining  a  nozzie  profile  is  the  calculac 
tion  of  the  right  side  of  equation  (10).    All  of  the  values  used  are  those  ob- 
tained from  equilibrium  flow  with  the  exception  of  the  rate  constants.     The 
rate  constants,  however,   may  be  assumed  to  be  functions  of  the  equilibrium 
temperature  at  each  point  in  the  nozzle.     Since  equilibrium  flow  is  indepen- 
dent of  the  nozzle  geometry,  the  same  values  calculated  for  one  system 
(chamber  conditions,  O/F  ratio,  and  reactants)  can  be  used  in  the  study  of 
any  nozzle  geometry,  regardless  of  size. 

Values  for  d  Qn0Cjj/fn)  ]  /d&    may  be  obtained  graphically  or  by  analyt- 
ical   methods  onfje  Xjj  and  m  are  known  for  different  area  ratios.    This  per- 
mits the  right  side  of  equation  (10)  to  be  calculated  and  plotted  versus  freeze 
area  ratio,  and  now  represents  B/(d6/dL)  as  a  function  of  freeze  area  ratio. 

To  convert  this  into  an  axisymmetric  nozzle  profile,  the  following 

geometric  relationship  may  be  used  [~7]  : 

d6     m    d(r/r*)2    =    2r_    d(r/r*)  (ll) 

dL  dL  '  r*  dL 

and  dr/dL  =  tan  0  (12) 

where  9  is  the  half -angle  of  the  nozzle.     Then 

dG      .    2r   tan  8    .     2  V6^tan  9  (13) 

dL  r*       r*  r* 

or  tan  8    =    <j6   -jg-  (14) 

dL    2f6^ 

Since  d£/dL  is  available  from  the  previous  calculation  of  B/(de  /dL), 
equation  (14)  provides  local  values  of  0  for  every  freeze  area  ratio.    It  should 
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be  noted  that  nozzle  size  is  now  present-  in  that  equation  (14)  involves  the 
throat  radius,   r*. 

The  half -angle,   9,   is  conveniently  calculated  by  obtaining  an  expres- 
sion for  B/(de/dL)  as  a  function  of  area  ratio.     If  B  is  assumed  to  be  unity, 
then  a  polynomial  form  for  d£/dL  is 

dL/de   =    bb  +  bx€   +  b^2  +  b^3  +  b4£4  (15) 

where  b0,   b,,   etc.  are  appropriate  coefficients. 

Equation  (15)  may  be  integrated  with  the  boundary  conditions  that 

e  =    1.0  at  L   =   O. 

2  3  4  c 

L   =    boe    +    —i —  *    -2 +    — ^ +  -4 (16) 

°  2  3  4  5 

-  (bQ  +  bx  +  b2  +  b3  +  b4) 

Using  equations  (14),   (15),  and  (16),  the  half -angle  and  nozzle  length 

may  be  calculated  for  each  freeze  area  ratio.    Since  the  throat  radius  has 

been  set,  the  throat  area  is  known.    With  the  area  ratio  and  the  throat  area, 

the  local  area  and  thus  the  nozzle  radius  are  determined.     If  an  equilibrium 

margin  is  desired,  it  must  be  applied  to  both  equations  (15)  and  (16). 

Special  provisions  must  be  made  near  the  nozzle  throat.    The  right  side 

of  equation  (10)  becomes  very  small  as  the  area  ratio  approaches  one.    Thus, 
d£/dL  is  large  and,  as  equation  (14)  shows,   the  half -angle  approaches  90 
degrees.     This  is  unacceptable  from  a  gas  dynamics  approach.    White  this 
difficulty  is  not  the  concern  of  this  study,  it  should  be  recognized  as  an 
area  requiring  special  handling. 

Some  method  of  fairing  in  the  profile  near  the  throat  must  be  used  until 
the  area  ratio  permits  the  contour  determined  by  equation  (14).    When  this 
point  in  the  nozzle  is  reached,  the  values  of  L  must  be  adjusted  since  they 
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no  longer  conform  to  the  calculated  contour.     This  correction  is  a  constant 
shift  of  the  L  values.    As  long  as  the  faired  section  near  the  throat  has  local 
half-angles  less  than  that  given  by  equation  (14),  the  flow  should  remain  in 
near-equilibrium.     Since  the  allowable  half-angles  in  this  region  are  large, 
no  particular  difficulty  is  encountered. 

At  large  area  ratios,   0  becomes  very  small.     Thus,   it  is  usually  desir- 
able to  terminate  the  contouring  after  the  majority  of  the  benefits  of  increased 
recombinationihave  been  realized.     Further  expansion  may  be  accomplished 
at  more  rapid  rates  with  consequent  near-frozen  flow. 
4.         Results  and  Discussion. 

A  study  of  a  particular  system  was  performed  to  check  the  validity  of 
the  proposed  method  of  nozzle  design.    The  hydrogen -oxygen  system  was 
chosen  because  of  its  relative  simplicity  and  limited  number  of  combustion 
products.    A  chamber  pressure  of  300  psia  and  an  oxidizer-to-fuel  weight 
ratio  of  six  was  chosen  as  representative  of  current  combustion  conditions. 

Equilibrium  and  frozen  values  for  the  above  mentioned  conditions  were 
computed  with  the  aid  of  a  computer  program  written  by  Gordon  and  Zeiznik 
[jj,   63  .     Tabulated  results  may  be  found  in  Appendix  II.     The  computer  pro- 
gram was  written  for  the  IBM  7090  system.     Since  the  computer  at  the  Naval 
Postgraduate  School  is  a  CDC  1604,   several  modifications  due  to  system 
differences  were  necessary  so  that  the  equilibrium  calculations  could  be 
made.     The  basic  program  was  unchanged. 

A  least-squares  curve  fitting  routine  was  used  to  obtain  a  third  order 
polynomial  relating  In  (Xjj/m)  and  area  ratio. 

23 


In  (XH/ffi)    =   aQ  +  axS  +  a2^2  +  a363  (17) 

where:         aQ    =    5.22823 

a-    =    5.22515  x  10"1 

a2    =    -4.27257  x  10"2 

a3    =    1.61956  x  10"3 
Thend  [in  (XH/!H)]  /d&    =    ax  +  2a2£  +  3a3e2  (18) 

The  rate  constants  k^,  k3,  and  k.4  were  calculated  from  the  following 
relationship: 

k    =    DTE  exp  (F/T)  (19) 

where  T  is  the  local  equilibrium  temperature  and  D,   E,   and  F  are  constants 

listed  in  Table  I  pL8]. 

TABLE    I 
REACTION  RATE  CONSTANTS 


Reaction 

D 

E 

F 

1 

2.3555  x  1017 

=1.4998 

-1493 

3 

5.5000  x  1011 

0.0000 

-10000 

4 

1.7268  x  1016 

-1.0900 

5375 

Using  equation  (18),   equation  (19),  and  the  equilibrium  values,  the 
right  side  of  equation  (10)  was  then  calculated.    The  results  are  given  in 
Table  II.    As  seen  from  equation  (10),  these  values  also  represent  B/(d€/dL). 
Thus,  the  locus  of  points  shown  in  Fig.   2  (EM  =1.0)  defines  the  hypothesized 
boundary  below  which  departures  from  near-equilibrium  flow  probably  occur. 
The  values  for  EM  =1.5  were  obtained  by  multiplying  the  basic  curve 
(EM    =    1.  0)  by  1.  5.     The  choice  of  1.  5  was  arbitrary.     In  view  of  the 
possible  inaccuracies  introduced  by  the  simplifying  assumptions  made  in 
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the  formulation  of  the  problem,   it  was  assumed  a  factor  of  approximately 
EM  =1.5  was  sufficient.     It  should  be  remembered  that  EM  =1.0  represents 
the  point  where,  theoretically,  the  flow  is  just  able  to  maintain  near-equili- 
brium.   A  large  value  for  EM  would  result  in  a  long  nozzle  since  the  expan- 
sion would  be  slower.    Also,  an  additional  important  aspect  of  this  design 
method  is  to  determine  the  maximum  rate  of  expansion  which  may  be  used 
and  still  maintain  near-equilibrium  conditions. 

A  polynomial  of  B/(d£/dL)  as  a  function  of  area  ratio  was  obtained 
with  a  curve-fit  routine  which  resulted  in  the  following  coefficients: 


bo 

=    9.28458xl0~3 

bl 

=  -2.39544  xlO"2 

b2 

=    1.34259  xlO"2 

b3 

=    3.95404xl0'3 

b4    =-3.83076  xlO"5 

Equation  (14)  was  used  to  determine  the  maximum  expansion  angle 
that  the  flow  could  tolerate  without  departing  from  near-equilibrium.     These 
values  are  shown  in  Table  III.     It  should  be  noted  that  small  area  ratios  yield 
large  values  for  0„ 

Calculations  were  performed  with  nozzle  throat  radii  of  2.  0  inches, 
0.96  inches,  and  0.  12  inches,,    These  different  size  nozzles  were  identified 
as  Run  1,   Run  2,  and  Run  3.    A  designation  was  added  to  each  run  to  indicate 
whether  the  nozzle  was  a  straight  conical  expansion  (S),  a  contoured  expan- 
sion (C),  or  a  contour  expansion  with  EM  =1.5  (CM).     For  example,  the  run 
with  the  two-inch  throat  that  had  a  contoured  expansion  with  EM  =1.0  was 
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designated  Run  1-C. 

A  maximum  of  45  degrees  was  used  for  all  runs  except  Run  3 -CM  where 
8  was  limited  to  30  degrees.     In  all  cases,  a  circular  arc  of  radius  0.  4  times 
the  throat  radius  was  used  to  fair  in  the  profile  near  the  throat  p.3].    A  lin- 
ear section  was  used  on  the  contoured  nozzles  (1-C,    l=CMt   2-C,  and 
3 -CM)  which  was  carried  to  the  point  where  the  area  ratio  indicated  that  the 
contoured  expansion  angle  as  shown  in  Table  III  was  permissable.    The 
straight  conical  expansions  were  all  at  15  degrees .    See  Table  IV  for  a  fur- 
ther description  of  the  nozzles.     Figs.  3,   4,   5,  and  6  show  the  profiles  used. 
Tables  V  to  XI  give  the  nozzle  geometry  for  each  run. 

Exact  calculations  for  all  seven  runs  were  made  with  the  computer 
program  written  by  Zupnik,   Nilson,  and  Sarli  [18].    As  with  the  equilibrium 
computer  program,   certain  system  modifications  were  necessary  to  convert 
from  the  IBM  7090  to  the  CDC  1604.    Tabulated  values  obtained  in  these  cal- 
culations are  shown  in  Appendix  III. 

The  benefits  of  maintaining  equilibrium  flow  can  be  seen  from  Fig.  7. 
Even  at  the  throat,  the  specific  impulse  difference  is  five  sees.    As  the  ex- 
pansion continues,  the  difference  increases  to  20  sees  at  an  area  ratio  of  40. 

Figs.  8,  9,  and  10  illustrate  the  effect  of  contouring  on  performance 
for  the  different  nozzle  sizes.    The  contour  is  most  effective  on  run  3 -CM 
and  is  insignificant  on  run  2-C.     The  apparent  explanation  for  run  2-=C  is 
that  this  run  is  not  modified.     That  is, '  it  does  not  have  the  equilibrium  margin 
that  runs  1-CM  and  3-CM  have.    Also,  a  small  nozzle  is  subject  to  earlier 
freezing  than  is  a  large  one.    This  is  due  to  the  smaller  residence  time  in  a 
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small  nozzle  which  gives  less  time  .tor  the  recombination  reactions  to  take 
place. 

Runs  2-C  and  2-S  illustrate  the  effect  of  any  contour  on  nozzle  length. 
Fig.   12  shows  that  a  specific  impulse  of  365  lb-sec/lb  may  be  obtained  with 
a  0.  10  foot  nozzle  for  the  contoured  nozzle.     The  straight  conical  nozzle, 
though,  needs  0.  17  feet  for  the  same  performance.     For  a  given  length  of 
0o08  feet,  the  specific  impulse  difference  is  28  sees. 

The  most  significant  result  is  that  shown  by  the  small  nozzle.     Even 
at  relatively  small  area  ratios  the  effect  of  contouring  is  evident.    At  an  area 

ratio  of  40,   the  difference  in  specific  impulse  is  six  sees.    While  this  is  not 

■ 

a  large  percentage  increase,  the  effect  over  a  long  period  would  be  appreci- 
able,,   Thus,   long  burn  times  would  require  less  fuel  for  the  same  total  per- 
formance.   Where  weight  is  a  critical  factor,   contouring  is  particularly  ad- 
vantageous.   For  a  hotter  system  such  as  hydrogen^fluorine,  the  percentage 
increase  in  performance  is  expected  to  be  more  significant. 

The  variation  of  hydrogen  concentration  is  given  in  Figs.    13,    14,  and 
15.    The  equilibrium  value  is  seen  to  decrease  uniformly  with  area  ratio. 
In  each  nozzle  size,  though,  the  exact  concentrations  show  a  leveling  off 
which  indicates  a  departure  from  equilibrium.     For  the  conical  nozzles  the 
change  is  gradual,   but  it  occurs  where  the  allowable  expansion  angle  is 
exceeded  (15  degrees).    All  of  the  contoured  nozzles  shown  a  marked  "knee" 
in  the  hydrogen  concentration  when  the  contouring  is  stopped  (a bout  can  areata tio 
of  4.  0).     Thus,   the  effect  of  contouring  a  nozzle  to  delay  departure  from 
equilibrium  appears  to  be  significant. 
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Each  of  the  figures  illustrating  hydrogen  concentration  shoves. a, gradual 
deviation  from  the  equilibrium  values  at  small  area  ratios.    It  is  felt  that  this 
is  the  effect  of  finite  rate  constants  and  not  due  to  inaccuracies  in  the  exact 
calculations.    When  large  rate  constants  (k  =  lO2^)  were  used,  the  exact 
values  essentially  duplicated  those  of  the  equilibrium  calculations.    Also, 
when  small  rate  constants  were  used  (k  =  100),  the  eoncentratttans  were 
invariant  indicating  frozen  or  non-reacting  flow. 

The  temperature  distribution  during  expansion  is  given  in  Fig.   16. 
Run  3-CM  has  larger  temperatures  than  Run  3-S.    This  indicates  that  more 
recombination  is  occurring  in  the  contoured  nozzle  than  in  the  straight  con- 
ical nozzle. 

As  mentioned  above,  the  amount  of  equilibrium  margin  used  was  arbi- 
trarily chosen  as  1.5.    Further  investigation  should  indicate  an  optimum 
factor  for  EM.    A  further  improvement  would  be  to  correlate  the  increase  in 
weight  due  to  a  longer  nozzle  with  the  increased  performance. 

The  approach  used  in  this  thesis  was  formulated  with  hydrogen  dis- 
appearance as  the  indication  of  the  freeze  point  and  subsequent  deviation  from 
near-equilibrium  conditions.    Further  study  might  show  that  other  species 
such  as  OH  or  monatomic  oxygen  give  better  results.     Possibly  a  combina- 
tion of  species  could  be  used.    Sarli,   et  al.  p.6]  investigated  this  modified 
approach  in  determining  the  freeze  point  for  multiple  reactions.    This  could 
serve  as  a  significant  refinement  to  the  basic  method  of  determining  optimum 
nozzle  contours. 
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5.        Conclusions. 

Optimum  contouring  of  a  nozzle  by  use  of  a  sudden  freeze  criterion 
results  in  a  delay  of  the  departure  from  near-requilibrium  flow  and  an  increase 
in  performance.    Xhis  is  relatively  more  effective  with  smaller  nozzles  where 
recombination  normally  lags  because  of  the  short  residence  time. 

Direct  application  of  the  Bray  criterion  for  sudden  freezing  does  not 
yield  the  best  results;  however,  an  equilibrium  margin  applied  to  the  basic 
curve  increases  the  performance.    Further  study  in  this  area  should  indicate 
an  optimum  equilibrium  margin. 

Hydrogen  may  not  be  the  ideal  species  to  indicate  the  freeze  point 
in  determining  the  optimum  nozzle  profile  for  this  system.    Other  possibilities, 
including  combinations  of  species ,   should  be  investigated. 
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TABLE  II 
FREEZE  POINT  LOCATION 


Area 

B 

Area 
Ratio 

B 

Ratio 

d    /dL 

d    /dL 

1.083 

. 0036307 

1.903 

.039031 

lo  124 

.0045955 

2.079 

.052248 

1.168 

.0056646 

2.251 

.067470 

1.213 

.0068524 

2.420 

.084753 

lo259 

.0081592 

3.976 

.365240 

1.353 

.0111240 

5.902 

1. 102000 

1,539 

.0184950 

7.891 

2.687200 

lo723 

.0277860 

13.355 

39.807000 

TABLE  III 
MAXIMUM  EXPANSION  ANGLE 
Half-angle  (Degrees) 


Area 

Run 

Run 

Run 

Run 

Ratio 

1-C 

1-CM 

2-C 

3-CM 

U04 

87.69 

86.54 

85.20 

44.76 

1.  15 

86.00 

84.01 

81o70 

29.75 

lo40 

79.73 

74.79 

69.31 

12.44 

U77 

64.03 

53.86 

44.59 

4.70 

U92 

56.21 

44o89 

33.65 

3.42 

2ol4 

44.75 

33.46 

25.45 

2.27 

2.21 

41.35 

30.40 

22.90 

2.02 

2o35 

35.12 

25.  12 

18.66 

1.61 

2.50 

29.38 

20„58 

15.13 

1.29 

2.80 

20.68 

14.13 

10.27 

0.87 

3,20 

13.39 

9.02 

6.52 

0.55 

4.00 

6.40 

4„28 

3.08 

0.26 

5,00 

3.07 

2.05 

1.47 

0.12 
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TABLE  IV 
NOZZLE  PARAMETERS 


Throat 

Throat 

Radius 

Area 

Run 

(ft)  _ 

.Jit2) 

l-S  0.  167  0.  087  3  Circular  arc  to  15°  straight  expan- 

sion. 

1-C  0o  167  0.  0873  Circular  arc  to  45°  ;  then  45°  straight 

to  45°  allowable  contour;  then  con- 
toured with  EM  =  1 .  0  to  point  of 
coincidence  with  Run  l-S;  then  15° 
straight  expansion* 

1-CM  0.  167  0.  087  3  Same  as  Run  1-C  except  contour  is 

with  EM  =  1.5. 

2-S  0o  08  0.0201  Same  as  Run  l-S  except  different 

size  throato 

2-C  0.08  0.0201  Same  as  Run  1-C  except  different 

size  throat. 

3-S  0.01  0.003142  Same  as  Run  l-S  except  different 

size  throat. 

3-CM  0. 01  0. 003142         Circular  arc  to  30°  ;  then  30°  straight 

to  30°  allowable  contour;  then  con- 
toured with  EM  =  1 .  5  to  area  ratio 
of  4.0;  then  15°  straight  expansion. 
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TABLE  V 
NOZZLE  GEOMETRY  FOR  RUN  1-S 


Lenqth 
(ft) 

Radius 

Area 
(ft2) 

Area 
Ratio 

0.000 

.0..  167 

0.,087  3 

1.000 

0.204 

0.219 

0. 1504 

1.723 

0.230 

0.226 

0. 1601 

1.834 

0.282 

0.240 

0.1805 

2.068 

0.330 

0.253 

0.2005 

2.297 

0.399 

0.271 

0.2308 

2.654 

0.461 

0.288 

0.2600 

2.979 

0.506 

0.300 

0.2825 

3.237 

0.577 

0.319 

0.3196 

3.662 

0.629 

0.333 

0.3484 

3o992 

0,697 

0.351 

0o3879 

4,444 

0.817 

0.384 

0.4631 

5.306 

1.103 

0.461 

0.668 

7.762 

1.294 

0.513 

0.8254 

9.458 

1.747 

0.634 

1.2634 

14,476 

2.160 

0.745 

1.7438 

19.982 

2.553 

0.851 

2.2734 

26.050 

3.261 

1.042 

3.4081 

39.052 
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TABLE  VI 
NOZZLE  GEOMETRY  FOR  RUN  1-C 


Length 
(ft) 

Radius 
(ft) 

Area 
(ft2) 

Area 
Ratio 

0.000 

0.167 

0,0873 

loOOO 

0.081 

0.219 

0.1504 

1,723 

0.101 

0.234 

0.1715 

U965 

0.110 

0.249 

0,1950 

2.234 

0.139 

0.267 

0.2234 

2.560 

0.187 

0.283 

0.2514 

2.881 

0.272 

0.301 

0.2839 

3.253 

0.328 

0.309 

0.3003 

3.441 

0.428 

0.322 

0.3251 

3.726 

0.553 

0.334 

0.3511 

4.023 

0.685 

0.354 

0.3943 

4.518 

0.840 

0.395 

0.4896 

5.610 

1.002 

0.439 

0.6052 

6.934 

1.199 

0.492 

0.7606 

8,715 

1.563 

0.589 

lo0907 

12.497 

2.276 

0.780 

1.9132 

21.923 

2o727 

0.901 

2.5514 

29.235 

3.692 

1.160 

4.2247 

48.410 
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TABLE  VII 
NOZZLE  GEOMETRY  FOR  RUN  1-CM 


Length 

Radius 

Area 

Area 

(ft)     , 

(ftL_ 

0.167 

(ft2) 
0.0873 

Ratio 

0.000 

1.000 

0o080 

0.219 

0.1507 

1.727 

0.093 

0.232 

0. 1694 

1.941 

0,102 

0.240 

0.1804 

2.106 

0.  115 

0.248 

0.1936 

2.210 

0.137 

0.259 

0.2108 

2,414 

0.  182 

0*274 

0.2366 

2.711 

0.312 

0.301 

0.2838 

3.251 

0.449 

0.318 

0.3169 

3.631 

0.640 

0.336 

0.3539 

4.055 

0,701 

0.352 

0.3897 

4.466 

0.837 

0.388 

0.4740 

5o431 

0.961 

0.423 

0.5609 

60  427 

1,263 

0.505 

0.8005 

9.173 

1.836 

0.658 

1.3610 

15. 594 

y 

2.699 

0.890 

2.4886 

28.516 

3.591 

1.  131 

4.0184 

46.046 
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TABLE  VIII 
NOZZLE  GEOMETRY  FOR  RUN  2-S 


Length 
(ft) 

Radius 
(ft) 

Area 
(ft2) 

Area 
Ratio 

0.000 

0.080 

0.0201 

1.000 

0o023 

0.085 

0.0226 

1.124 

0.035 

0.088 

0.0245 

1.217 

0.069 

0.097 

0.0298 

1.481 

0.103 

0.107 

0.0356 

1.772 

0.149 

0.119 

0.0442 

2.198 

0,170 

0.124 

0.0486 

2.417 

0.198 

0.132 

0.0546 

2.717 

0.246 

0.145 

0o0658 

3.270 

0.275 

0.153 

.    0.0731 

3.636 

0.312 

0.163 

0.0826 

4.125 

0.349 

0.172 

0.0932 

4.635 

0.422 

0.192 

0.1155 

5.746 

0.505 

0.214 

0. 1440 

7.160 

0.673 

0.259 

0.2108 

10.482 

0.916 

0.324 

0.3303 

16.422 

1.255 

0.415 

0.5412 

26.911 

1.463 

0.471 

0.6963 

34.625 

1,631 

0.516 

0.8356 

41.552 
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TABLE  IX 
NOZZLE  GEOMETRY  FOR  RUN  2-C 


Length 
(ft) 

Radius 
(ft) 

Area 
(ft2) 

Area 
Ratio 

0.000 

0.080 

0.0201 

1,000 

0.017 

0.085 

0.0226 

lol23 

0.024 

0.091 

0.0260 

1.295 

0.030 

0.097 

0.0296 

1.473 

0.038 

0.  105 

0o0346 

1.722 

0.051 

0.  114 

0.0414 

2,046 

0.067 

0.  122 

0o0465 

2.314 

0.095 

0.130 

0.0527 

2,621 

0,156 

0.140 

0.0616 

3,062 

0.234 

0.149 

0.0693 

3.445 

0.308 

0.161 

0.0819 

4.073 

0.356 

0.175 

0.0958 

4.764 

0.442 

0.198 

0.  1225 

6.093 

0.624 

0.246 

0. 1895 

9.425 

0.801 

0.293 

0.2700 

13,424 

0.981 

0.342 

0.3667 

18.236 

1.320 

0.432 

0.5875 

29,216 

1.528 

0.488 

0.7486 

37.160 

1.669 

0.526 

0.8689 

43,092 
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TABLE  X 
NOZZLE  GEOMETRY  FOR  RUN  3-S 


Length 
x  102 

(ft) 

Radius 
x  102 

(ft) 

1.000 

Area 
x  103 
(ft2) 

0.3142 

Area 
Ratio 

0.000 

lo000 

0.204 

1.041 

0,3402 

lo083 

0.432 

1.101 

0.3810 

1.212 

0.808 

1.202 

0.4538 

1.444 

1,607 

1.417 

0.6308 

2.008 

2.368 

1.620 

0.8246 

2.625 

2,992 

1.787 

1.0035 

3.194 

3.935 

2,040 

1.3079 

4,163 

5.991 

2.591 

2.1093 

6,714 

9.933 

3.647 

4.1794 

13.303 

14.211 

4.794 

7.2194 

22.980 

18.490 

5o940 

11.0843 

35,282 
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TABLE  XI 
NOZZLE  GEOMETRY  FOR  RUN  3-CM 


Length 
x  102 

Radius 

x  102 

(ft) 

Area 

x  103 

(ft2) 

Area 
Ratio 

0.000 

1.000 

0.3142 

1.000 

0.176 

1.041 

0.3402 

1.090 

0.371 

1.132 

0.4023 

1.281 

0.563 

1.183 

0.4400 

1.401 

1.772 

1.334 

0.5593 

1.780 

3.199 

1.423 

0.6358 

2.024 

6.940 

1.553 

0.7576 

2.412 

12.555 

1.668 

0.8738 

2.781 

29,700 

1.810 

1.0289 

3.275 

39.992 

1.917 

1.1547 

3.676 

52,417 

2.006 

1.2639 

4.023 

57.021 

2.147 

1.4488 

4.612 

57.671 

2.325 

1.6981 

5.406 

59.936 

2.985 

2.7995 

8.911 

63.627 

3.961 

4.9299 

15.692 

66.395 

4.705 

6.9552 

22. 139 

68.324 

5.222 

8.5669 

27.270 

70.  170 

5.717 

10.2664 

32.679 

72.015 

6.211 

12.1190 

38.576 

72.938 

6.458 

13. 1030 

41.708 
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FIGURE  2 


NOZZLE  GEOMETRY  PARAMETERS  AS  A 
/  /  FUNCTION  OF  FREEZE  AREA  RATIO 
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FIGURE  10 
SPECIFIC  IMPULSE 
COMPARISON  OVER  ENTIRE 

NOZZLE  FOR  THROAT 
RADIUS  OF  0.  12  INCHES 
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FIGURE  16 
TEMPERATURE  DISTRIBUTION 
FOR  THROAT  RADIUS  OF  0.  12  INCHES 
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APPENDIX    I 
EQUATIONS  FOR  EXACT  CALCULATIONS 

The  equations  involved  in  the  analysis  of  nozzle  flow  with  a  chemical 
reaction  include  those  of  gas  dynamics  and  of  reaction  kinetics  [181  . 

Gas  Dynamic  Equations 

Overall  mass  conservation: 

ffynr   =  rh  (1-1) 

Total  energy  conservation: 

h  hin  +~£   =  Ht>  (1-2) 

where: 


hi  '  £Cec  dT  +  (Hu^ 


Momentum  conservation: 

nrjor  +  fy  dp  -  0  (M) 

Equation  of  state: 

P^/OfiTZfT  (1=4) 

Speed  of  sound: 
where    /  is  defined  as: 


y  ==  £l£*ig d-7) 


Mach  number: 


M   =    ^  ft. 
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Reaction  Kinetics  Equations 
General  form  of  nr  chemical  reactions: 

z.  '^V  '  *     ^    .  "V  (1=9) 

where  24/    and  3B$£     are  the  stoichiometric  coefficients  of  the  ft1  species 
in  the  jth  reaction,  reactants  and  products  respectively,  and  k.£j  and  k^j 
are  the  specific  reaction  rates  of  the  forward  and  backward  chemical  changes 
in  the  fi1  reaction  of  the  chemical  system. 

The  rates  of  change  of  the  concentration  of  the  species  M^  are  deter- 
mined by  the  relation: 


The  atomic  species  continuity  equations  are: 


(+z  IjZ,  •'";  n?) 


U-io) 


(i- 11) 


where  B^  is  the  number  of  the  i*"  atomic  species  in  the  j1*1  molecular  species, 
n  designates  the  total  number  of  chemical  species,  and  na  is  the  number  of 
atomic  species. 

The  molecular  species  continuity  equations  are: 


M  C*i  •**)  fa  g;  ("V* 

7*  i  W  **l 


""  A/  $*  (^**)     '  J l  r  n>H> r»+s> . . . 


a-12) 
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It  should  be  noted  that  of  the  n  equations  represented  by  (1=11)  and 
(1=12),  only  n-1  are  independent  because  of  the  overall  mass  conservation 
equation  (1=1).    The  exact  calculations  involve  the  simultaneous  solution 
of  equations  (1=1),   (1=2),    (1=3),   (1=4),   (I- 10),   (1=11),  and  (1-12). 
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APPENDIX    II 

EQUILIBRIUM  AND  FROZEN  RESULTS 

TABLE    n-1 
Equilibrium  Composition  Values 


Area 
Ratio 

Pressure 
(Atmos) 

20.414 

Pressure 
Ratio 
Oc/p) 

1.000 

Temp           ( 
Deg  K 

3439           ] 

lamina 
L.  1313 

Specific 

Heat 

Cal/gm-°K 

2.9563 

Enthalpy 
Cal/gm 

Chamber 

0 

Throat 

11.811 

1.728 

3274           ] 

L.  1294 

2.7698 

-277.8 

1.083 

8.165 

2.500 

3166           ] 

L.1288 

2.6193 

-455.7 

1.124 

7.423 

2.750 

3139           ] 

L.  1288 

2.5773 

-500.4 

1.  168 

6.805 

3.000 

3114            ] 

..1288 

2.5380 

-540.8 

1.213 

6.281 

3.250 

3091            ] 

L.1288 

2.5010 

-577.6 

1.259 

5.832 

3.500 

3070           ] 

L.  1289 

2.4660 

-611.3 

1.353 

5.103 

4.000 

3033           ] 

l.  1291 

2.4013 

-671.5 

1.539 

4.083 

5.000 

2971            ] 

L.  1297 

2.2887 

-770.0 

1.723 

3.402 

6.000 

2920           ] 

..1304 

2. 1932 

-848.8 

1.903 

2.916 

7.000 

2877            ] 

L.  1313 

2.1105 

-913.8 

2.0Z9 

2.552 

8.000 

2840           ] 

L.  1322 

2.0376 

-969.3 

2.251 

2.268 

9.000 

2807           ] 

l.  1332 

1.9727 

-1017.6 

2.420 

2.041 

10.000 

2778           ] 

l.1342 

1.9144 

-1060.2 

3.976 

1.021 

20.000 

2582           ] 

L.1442 

1.5364 

-1327.5 

5.902 

0.600 

34.020 

2425            ] 

L.1567 

lo2786 

-1517.3 

7.891 

0.408 

50.000 

2306            ] 

L.  1682 

1.1237 

-1646.5 

13.355 

0.204 

100.000 

2083           ] 

L.  1913 

0.9259 

-1861.6 

45.521 

0.041 

500.000 

1573            ] 

L.  2311 

0.7515 

-2275.2 

77.299 

0.020 

5000.000 

989           ] 

L.2437 

0.7190 

-2419.0 

sa 


Area 
Ratio 


TABLE  II -2 

Equilibrium  Composition  Values 

Mean  Mach  C*  C„ 

Mol.  wt.      Number       Ft/sec 


Specific  Impulse 
(lb- sec/lb) 

1sp)v  Xsp 


Chamber     13.043        0.000 


Throat 

13.223 

1.000 

7655 

.654 

293.  1 

155.5 

1.083 

13.341 

1.308 

7655 

.837 

302.2 

199.1 

1.124 

13.370 

1.379 

7655 

.877 

305.9 

208.7 

1.  168 

13.397 

1.440 

7655 

.912 

309.6 

216.9 

1,213 

13.421 

1.495 

7655 

.942 

313.0 

224.2 

1.259 

13.444 

1.545 

7655 

.969 

316.3 

230.6 

1.353 

13.484 

1.631 

7655 

1.016 

322.2 

241.7 

1.539 

1)0.548 

1.769 

7655 

1.088 

332.1 

258.8 

1,723 

13.599 

1.876 

7655 

1.142 

340.0 

271.7 

1.903 

13.641 

1.963 

7655 

1.  185 

346.6 

282.0 

2.079 

13.677 

2.037 

7655 

1.221 

352.2 

290.4 

2,251 

13.707 

2.101 

7655 

1.251 

357.  1 

297.6 

2.420 

13.733 

2.  157 

7655 

1.277 

361.3 

303.7 

3.976 

13.887 

2.506 

7655 

1.429 

387.2 

339.9 

5,9q2 

13.978 

2.759 

7655 

1.527 

404.6 

363.3 

7.891 

14.028 

2.938 

7655 

1.591 

416.  1 

378.5 

13.355 

14.082 

3.261 

7655 

1.692 

434.2 

402.5 

45.521 

14. Ill 

4.085 

7655 

1.870 

466.6 

444.9 

77.299 

14.112 

4.479 

7655 

1.928 

477.2 

458.8 

60, 


TABLE  II-3 


Frozen  Composition  Values 


Area 
Ratio 

Pressure 
(Atmos) 

20.41:4 

Pressure 
Ratio 

(Po/P) 
1.000 

Temp 
Deg  K 

Gamma 
1.2057 

Specific 

'Heat 
Cal/qm-°K 

.8932 

Enthalpy 
Cal/gm 

Chamber 

3439 

0 

Throat 

11.483 

1.778 

3115 

1.2096 

.8792 

-287.1 

1.065 

8.165 

2.500 

2936 

1.2123 

.8702 

-444.2 

1.101 

7.423 

2.750 

2887 

1.21(30 

.8676 

-486.5 

1.139 

6.805 

3.000 

2843 

1.2137 

.8652 

-524.5 

1.179 

6.281 

3.250 

2803 

1.2144 

.8629 

-558.9 

1,219 

5.832 

3.500 

2767 

1.2151 

.8608 

-590.4 

1.301 

5.103 

4.000 

2702 

1.2162 

.8570 

-646.0 

1.464 

4.083 

5.000 

2597 

1.2183 

.8503 

-736.1 

1.624 

3.402 

6.000 

2513 

1.2201 

.8447 

-807.0 

1.779 

2.916 

7.000 

2444 

1.2216 

.8398 

-865.2 

1.931 

2.552 

8.000 

2385 

1.2230 

.8355 

-914.3 

2.078 

2.268 

9.000 

2334 

1.2243 

.8316 

-956.7 

2.222 

2.041 

10,000 

2290 

1.2255 

.8281 

-993.8 

3.523 

1.021 

20.000 

2012 

1.2339 

.8037 

-1220.7 

5.098 

.600 

34.023 

1817 

1.2413 

.7837 

-1375.5 

6.702 

.408 

50.000 

1684 

1.2472 

.7686 

-1478.2 

11.048 

.204 

100.000 

1464 

1.2590 

.Z406 
1 

-1644.2 

35.866 

.041 

500.00 

1035 

1.2924 

.6735 

-1948.5 

59.600 

.020 

1000.000 

882 

1.3074 

.6480 

-2049.5 
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TABLE  II-4 
Frozen  Composition  Values 

Area  Mean  Mach  C*  Gp-  Specific  Impulse 

Ratio  Mol.  wt.      Number       Ft/sec  (lb-sec/lb) 

sp)v  sp 


Chamber  13.043  0 

Throat  13.043  1.000 

1.065  13.043  1.280 

1.101  13.043  1.350 

1.139  13.043  1.412 

1.179  13.043  1.468 

1.219  13.043  1.518 

1.301  13.043  1.606 

1.464  13.043  1.748 

1.624  13.043  1.859 

1.779  13.043  1.950 

1.931  13.043  2.028 

2.078  13.043  2.096 

2.222  13.043  2.156 

3.523  13.043  2.541 

5.098  13.043  2.830 

6.702  13.043  3.039 

11,048  13.043  3.422 

35.866  13.043  4.373 

59.600  13.043  4.830 


7473 

.680 

288.7 

158.0 

7473 

.846 

295.6 

196.6 

7473 

.886 

298.8 

205.7 

7473 

.920 

301.8 

213.6 

7473 

.949 

304.8 

220.5 

7473 

.976 

307.6 

226.6 

7473 

1.021 

312.6 

237.1 

7473 

1.090 

321.1 

253.1 

7473 

1.  141 

327.9 

265.0 

7473 

1.181 

333.4 

274.4 

7473 

1.214 

338.  1 

282.1 

7473 

1.242 

342.2 

288.5 

7473 

1.266 

345.7 

294.1 

7473 

1.403 

366.8 

325.9 

7473 

1.489 

380.8 

346.0 

7473 

1.544 

389.8 

358.6 

7473 

1.628 

403.9 

378.2 

7473 

1.773 

428.4 

411.8 

7  473 

1.818 

436.  1 

422.3 
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TABLE  II- 5 
EQUILIBRIUM  FLOW  SPECIES  CONCENTRATIONS 
Mol  Fractions  xlO3 


Area 
Ratio 

Ol 

Hl 

°2 

H2 

OH 

HzO 

Chamber 

8.16 

56.88 

7.75 

252.15 

62.82 

612.24 

Throat 

6.06 

49.35 

6.14 

248.38 

52.17 

637.90 

1.083 

4.81 

44.33 

5.07 

246.06 

45.20 

654.53 

1.124 

4.51 

43.03 

4.80 

245.50 

43.44 

658.71 

1,168 

4.25 

41.86 

4.56 

245.00 

41.84 

662.49 

1.213 

4.01 

40.77 

4.34 

244.56 

40.38 

665.93 

1.259 

3.80 

39.78 

4.14 

244.16 

39.04 

669.08 

1.353 

3.43 

37.98 

3.79 

243.48 

36.66 

674.67 

1.539 

2.85 

34.99 

3.22 

242.44 

32.76 

683.73 

1.723 

2.43 

32.56 

2.79 

241.69 

29.67 

690.87 

1.903 

2.10 

30.52 

2.44 

241.13 

27.12 

696.69 

2.079 

1.83 

28.76 

2.15 

240.72 

24.98 

701.56 

2.251 

1.61 

27.22 

1,92 

240.40 

23.14 

705.72 

2.420 

1.43 

25.85 

1.72 

240. 15 

21.53 

709.32 

3.976 

0.55 

17.19 

0.69 

239.55 

12.16 

729.86 

5.902 

0.21 

11.26 

0.26 

240. 17 

6.74 

741.36 

7.891 

0.08 

7.63 

0.11 

240.99 

3.94 

747.24 

13.355 

0.01 

3.04 

0.01 

242.56 

1.12 

753.25 

45.521 



0.10 



243.94 

0.01 

755.95 

77.299 

____ 

0.01 

M  i 

243.99 

,1  n  i , 

755.99 

u 


APPENDIX  III 

EXACT  RESULTS 

TABLE  III-l 

EXACT  FLOW  SPECIES  CONCENTRATIONS  FOR  RUN  1-S 

Species  Concentration 
(lb-mol/1^  x  104 

Area                 O                   H                  02                 H2                OH              H20 
Ratio  


1.723  1.79  23,91  2.05  177.83  21.83  507.98 

1.834  1.72  23.66  2.02  177.50  21.19  508.75 

2.068  1.51  22.45  1.81  176.51  19.27  511.30 

2.297  1.33  21.37  1.63  175.72  17.68  513.43 

2.645  1.14  20.10  1.41  174.77  15.77  515.96 

2.979  1.00  19.21  1.26  174.07  14.35  517.82 

3.237  0.92  18.69  1.18  173.63  13.46  518.96 

3.662  0.82  18.07  1.06  173.03  12.27  520.48 

3.992  0.77  17.73  1.00  172.64  11.53  521.41 

4.444  0.71  17.41  0.93  172.18  10.69  .522.45 

5.306  0.63  17,10  0.84  171.47  9.48  523.91 

7.762  0.52  17.15  0.72  170.10  7,48  526.25 

9.458  0.47  17.45  0.75  169.33  6.52  527.35 

14.476  0.40  18.42  0.61  167.81  4.86  529.23 

19,982  0.35  19.29  0.57  166.70  3.77  530.45 

26.050  0.31  20.01  0.54  165.84  2.99  531.33 

39,052  0.24  21,03  0.51  164.70  1,98  532.47 
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TABLE  III-2 

EXACT  FLOW  SPECIES  CONCENTRATIONS  FOR  RUN  1-C 

Species  Concentration 
(lb-mol/lb)  x  104 


Area 
Ratio 

O 

H 

^2 

H2 

OH 

HzO 

1.723 

1.79 

23.91 

2.05 

177.83 

21.80 

507.98 

1.965 

1.71 

23.93 

2.03 

177.06 

20.57 

509.37 

2.234 

1.64 

24.00 

1.99 

176.21 

19.25 

510..84 

2.560 

1.44 

22.90 

1.81 

175.29 

17.41 

513.23 

2.881 

1.22 

21.35 

1.55 

174.46 

15.66 

515.71 

3.253 

1.00 

19.51 

1.28 

173.69 

13.84 

518.31 

3.4*i 

0.90 

18.66 

1.  16 

173.37 

13.03 

519.46 

3.726 

0.78 

17.49 

1.00 

172.97 

11.93 

520.99 

4.023 

0.67 

16.43 

0.87 

172.64 

10.93 

522.36 

4.518 

0.59 

15.79 

0.78 

172.18 

9.92 

523.64 

5.610 

0.52 

15.67 

0.70 

171.39 

8.64 

525.14 

6.934 

0.47 

15.77 

0.65 

170.64 

7.56 

526.37 

8.715 

0.43 

16.06 

0.60 

169.85 

6.55 

527.51 

12.497 

0.37 

16c79 

0.55 

168.63 

5.16 

529.07 

21.923 

0.29 

18.23 

0.49 

166.81 

3.34 

531.08 

48.410  0.19  20.09  0.44  164.74  1.49         533.15 
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TABLE  III-3 

EXACT  FLOW  SPECIES  CONCENTRATIONS  FOR  RUN  1-CM 

Species  Concentration 
(lb-mol/lb)  x  104 

Area                 O                  H                  02                H2                OH              H20 
Ratio  


3,251 

0.92 

18.70 

1.18 

3.631 

0.76 

17.15 

0,97 

1,727  1.78  23.80  2.05  177.94  21.95  507.86 

1.941  1.73  24,00  2.04  177.15  20.73  509.15 

2ol06  1.66  23.78  2.00  176.72  19.94  510.10 

2.210  1.57  23.30  1.96  176.21  18.98  511.32 

2.414  1.43  22o51  1.78  175.59  17.78  512.94 

2.711  1.23  21.16  1.55  174.77  16.09  515.28 

173.64  13.45  518.96 

173.07  11.92  521.07 

4,055  0.62  15.76  0.80  172.59  10.55  522.93 

4.466  0.59  15.68  0,77  172.21  9.92  523.66 

5,431  0.53  15„57  0.70  171.50  8.76  525.01 

6.427  0.48  15.63  0.66  170.90  7.88  526.01 

9.173  0.42  16.08  0.69  169.68  6.30  527.80 

15,594  0.34  17.27  0.52  167.90  4.37  529.95 

28,516  0.26  18.83  0.47  166.05  2.62  531.90 

46.046  0.19  19.91  0.44  164.87  1.57  533.06 


66 


TABLE  III-4 

EXACT  FLOW  SPECIES  CONCENTRATIONS  FOR  RUN  2-S 

Species  Concentration 
(lb-mol/lb)  x  104 

Area                O                  H                 02                H2               OH             HzO 
Ratio  


3.270 

1.15 

21.18 

1.48 

3.636 

1.08 

20.92 

1.40 

1.124  3.59  32.14  3.59  183.53  32o48  492.50 
1.217  3.10  30.93  3.43  182.45  30.69  495.08 
1,481  2.42  27.67  2.78  179.73  25.97  501.79 
1.772  1.98  25.39  2.34  177.83  22*54  506.55 
2.198  1.59  23o32  1.94  175.99  19.17  511.11 
2.417  1.45  22.63  1.80  175.30  17.90  512.78 
2.717  1.32  21.93  1.66  174.52  16.51  514.61 

173.42  14.60  517.04 

172.84  13.66  518.21 

4.125  1.00  20.76  1.32  172.19  12.64  519.46 
4.635  0.95  20.73  1.26  171.60  11.79  520.48 
5.746  0.87  20.91  1.18  170.57  10.40  522.12 
7.160  0.80  21.34  1.11  169.52  9.14  523.58 

10.482  0.71  22.48  1.04  167.74  7.20  525.76 

16.422  0.61  24.11  0.97  165.70  5.22  ,527.9% 

26.911  0.50  25.91  0.92  163.67  3.40  530.01 

34.625  0.44  26.73  0.90  162.79  2.64  530.87 

41.552  0.41  27.26  0.88  162.23  2.16  531.41 
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TABLE  III-5 

EXACT  FLOW  SPECIES  CONCENTRATIONS  FOR  RUN  2-C 

Species  Concentration 
>-mol/lb)  x  104 


Area                 O  H  02                H2  OH  H20 

Ratio  , , _____  _____     

1.123  3.59  32.14  3.59  183.53  32.48  492.50 

1.295  3,08  31.15  3.47  181.88  29.66  496.06 

1.473  2.79  30.05  3.23  180.39  27.14  499.35 

1.722  2.49  28.96  2.97  178.82  24.54  502.77 

2o046  2.18  27.69  2.67  177.25  21.94  506.28 

2.314  1.91  26.30  2.39  176.15  20.03  509.03 

2.621  1.63  24w63  2.07  175.11  18.11  511.86 

3.062  1.30  22o27  1.66  173.95  15.72  515.39 

3.445  1.07  20.46  1.38  173.20  14.00  517.92 

4.073  0.93  19.73  1.22  172.32  12.44  519.93 

4.764  0.87  19.79  1.15  171.57  11.37  521.20 

6.093  0.79  20.10  1.08  170.43  9.86  522.93 

9.425  0.69  21.23  0.98  168.44  7.60  525.49 

13.424  0.61  22.44  0.93  166.86  6.00  527.26 

18.236  0.55  23.56  0.90  165.54  4.76  528.64 

29.216  0.45  25.21  0.85  163.71  3.12  530.46 

37.160  0.40  25.96  0.84  162.92  2.43  531.23 

43.092  0.38  26.36  0.83  162.48  2.06  531.65 
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TABLE  III-6 

EXACT  FLOW  SPECIES  CONCENTRATIONS  FOR  RUN  3-S 

Species  Concentration 
(lbr-mol/lb)  x  104 


Area 

O 

H 

» 

°2 

H 

2 

OH 

H2 

O 

Ra 

tio 
083 

H 

1. 

3. 

80 

33. 

23 

3. 

80 

184. 

38 

33. 

87 

490. 

48 

1. 

212 

3. 

49 

33. 

24 

3. 

84 

183. 

11 

31. 

80 

492. 

78 

1. 

444 

3, 

17 

32. 

42 

3. 

66 

181. 

16 

28. 

46 

496. 

81 

2. 

008 

2. 

75 

31. 

94 

3c 

34 

178. 

03 

23. 

82 

502. 

49 

2. 

625 

2. 

55 

32. 

46 

3. 

19 

175. 

79 

20. 

89 

505. 

93 

3. 

194 

2. 

44 

33o 

22 

3. 

11 

174. 

22 

19. 

01 

508. 

07 

6. 

714 

2. 

15 

37. 

58 

2. 

99 

168. 

42 

12. 

85 

514. 

77 

13. 

303 

1. 

90 

42. 

12 

2. 

95 

163. 

47 

8. 

13 

519. 

81 

22. 

980 

1. 

72 

45. 

15 

2. 

95 

160. 

33 

5. 

24 

522. 

87 

35. 

282 

1. 

61 

46. 

84 

2. 

96 

158. 

61 

3. 

67 

524. 

53 
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TABLE  III-7 
EXACT  FLOW  SPECIES  CONCENTRATIONS  FOR  RUN  3-CM 


Species  Concentration 
(lb-mol/lb)  x  104 

Area 
Ratio 

O 

H 

o2 

H2 

OH 

H20 

1.090 

3.80 

33.23 

3.80 

184.38 

33.87 

490.48 

1.281 

3.49 

33.69 

3.86 

182.53 

31.00 

493.53 

1.401 

3.31 

33.08 

3.78 

181.59 

29.20 

495.67 

1.780 

2.61 

29.95 

3.16 

178.90 

24.56 

502.25 

2.024 

2.20 

27o79 

2.72 

177.48 

22.  16 

505.95 

2.412 

1.69 

24c  65 

2.  12 

175.76 

18.99 

510.82 

2.781 

1.34 

22.20 

1.70 

174.58 

16.57 

514.43 

3.275 

0.92 

18.53 

1.17 

173.41 

13.53 

518.96 

3.676 

0.79 

17.45 

1.01 

172.83 

12.18 

520.75 

4,023 

0.69 

16.54 

0.89 

172.42 

11.14 

522.14 

4.612 

0.66 

16.72 

0.85 

171.86 

10.39 

522.99 

5.406 

0.65 

17.33 

0.85 

171.18 

9.69 

523.70 

8.911 

0.62 

19.38 

0.85 

169.05 

7.53 

525.90 

15.692 

0.58 

21.61 

0.85 

166.77 

5.29 

528.18 

22. 139 

0.55 

22.81 

0.85 

165.55 

4.  12 

529.39 

27.270 

0.53 

23.44 

0.85 

164.91 

3.50 

530.02 

32.679 

0.51 

23.93 

0.85 

164.42 

3.04 

530.50 

38.576 

0.50 

24.33 

0.85 

164.02 

2.66 

530.89 

41.708 

0.50 

24.49 

0.85 

163.86 

2.50 

531.05 
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TABLE  III-8 


EXACT  FLOW  VALUES  FOR  RUN  1-S 


Area 
Ratio 

Velocity 
(ft/sec) 

Pressure 
(lb/ft2) 

Temp 
(Deg  R) 

Mach 
Number 

Specific      Density 
Impulse       x  103 
(lb-sec/lb)   (lb/ft3) 

340.1             12.06 

Gamma 
(Frozen) 

1.723 

8743 

7201 

5256 

1.82 

1.864 

1.834 

8964 

6497 

5180 

1.88 

344.2 

11.04 

1.156 

2o068 

9344 

5395 

5070 

1.98 

351o9 

9.39 

1.146 

2.297 

9645 

4613 

4981 

2.07 

358.1 

8.20 

1.148 

2.645 

10020 

3754 

4862 

2.17 

366,1 

6.85 

1.154 

2,979 

10313 

3164 

4760 

2.26 

372*4 

5.91 

1.159 

3.237 

10506 

2813 

4689 

2.32 

376.7 

5.34 

1.163 

3.662 

10780 

2364 

4581 

2.41 

382.7 

4.60 

1.169 

3.992 

10962 

2095 

4505 

2.48 

386.8 

4.15 

1.173 

4.444 

11178 

1805 

4411 

2.55 

391.6 

3.66 

1.178 

5.306 

11514 

1414 

4254 

2.68 

399.2 

2.97 

1.185 

7o762 

12134 

858 

3933 

2.93 

413.3 

1.96 

1.197 

9.458 

12455 

645 

3752 

3.08 

420.7 

1.54 

1.203 

14.476 

13027 

365 

3404 

3.37 

434.0 

0.96 

1.213 

19.982 

13407 

238 

3156 

3.60 

442.9 

0.68 

1.220 

26.050 

13690 

168 

2960 

3.78 

449.6 

0.51 

1.226 

39.052 

14078 

99 

2679 

4.08 

458.8 

0.33 

1.236 

Zl 


TABLE  III-9 


EXACT  FLOW  VALUES  FOR  RUN  1-C 


Area         Velocity  Pressure    Temp         Mach       Specific       Density   Gamma 
Ratio         (ft/sec)    (lb/ft2)     (Deg  R)     Mumber     Impulse           x  103     (Frozen) 
(lb-sec/lb)    (lb/ft3) 


1.723 


8743        7201 


5256 


1.82 


340.0 


12.05        1.577 


1.965 


9197         5803 


5088 


1.94 


348.7 


10.05        1.181 


2.234 


9589         4747 


4937 


2.06 


356.4 


8.48        1.163 


2.560 


9953        3886 


4819 


2.16 


364.  1 


7.13        1.141 


2.881        10238        3292 


4740 


2.25 


370.4 


6.16        1.129 


3.253        10512         2786 


4666 


2.33 


376.6 


5.31        1.126 


3.441        10633        2582 


4633 


2.37 


379.4 


4.96        1.127 


3.726        10800        2319 


4586 


2.42 


383.2 


4.51        1.131 


4.023        10956        2092 


4538 


2.47 


386.9 


4.12        1.144 


4.518        11189 


1783 


4442 


2.55 


392.1 


3.59        1.184 


5.610        11599         1321 


4242 


2.70 


401.3 


2.79        1.190 


6o934        11966 


988 


4051 


2.85 


409.6 


2.19        1.196 


8.715        12327 


725 


3850 


3.01 


417.9 


1.69        1.202 


12.497        12836 


21.923        13510 


29.235        13809 


48.410        14269 


447 

3547 

3.26 

429.7 

212 

3107 

3.65 

445.5 

145 

2898 

3.86 

452.6 

75 

2555 

4.23 

463.4 

1.13        1.210 


0.61        1.222 


0.45        1.229 


0.26        1.241 
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TABLE  III- 10 

EXACT  FLOW  VALUES  FOR  RUN  1-CM 

Area         Velocity  Pressure    Temp       Mach  Specific     Density   Gamma 

Ratio         (ft/sec)    (lb/ft2)     (Deg  R)    Number        Impulse      x  103       (Frozen) 

(lb-sec/lb)  (lb/ft3) 


1.727 

8752 

7169 

5252 

1.82 

340.2 

12.01 

1.201 

1.941 

9157 

5919 

5100 

1.93 

347.9 

10.22 

1.175 

2.106 

9354 

5363 

5032 

1.99 

351.8 

9.39 

1.164 

2.210 

9561 

4817 

4965 

2.05 

356.0 

8.56 

1.154 

2.414 

9793 

4251 

4893 

2.11 

360.9 

7.68 

1.145 

2.711 

10087 

3602 

4808 

2.20 

367o3 

6.64 

1.136 

3.251  10509  2798  4687  2.33  376.7  5.31  1.132 

3.631  10746  2408  4618  2.40  382.1  4.65  1.131 

4.055  10973  2074  4546  2,47  387.4  4.08  1.182 

4.466  11171  1810  4453  2.54  391.7  3,64  1.185 

5.431  11545  1380  4271  2.68  400.1  2.90  1.190 

6.427  11843  1096  4118  2.80  406.9  2.39  1.195 

9.173  12408  676  3805  3.05  419.8  1.60  1.203 

15.594  13122  333  3367  3.42  436.4  0.89  1.215 

28.516  13788  150  2915  3.84  452.1  0.46  1.228 

46.046  14231  80  2587  4.02  462.5  0.28  1.240 
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TABLE  III- 11 


EXACT  FLOW  VALUES  FOR  RUN  2-S 


Area 
Ratio 

Velocity 
(ft/sec) 

Pressure 
(lb/ft2) 

Temp 
(Deg  R) 

Mach 
Number 

Specific 

Impulse 

(lb-sec/lb) 

Density 

x  103 

(lb/ft3) 

Gamma 
(Frozen) 

1,124 

6715 

15707 

5650 

1.33 

306.0 

24.05 

1.734 

1.217 

7250 

13119 

5534 

1.46 

313.3 

20.57 

.1.134 

1.481 

8190 

9141 

5338 

1.68 

329.  1 

14.96 

1.138 

1.772 

8856 

6814 

5174 

1.85 

341.8 

11.57 

1.144 

2.198 

9525 

4891 

4980 

2.04 

355.3 

8.67 

1.154 

2o417 

9787 

4248 

4895 

2oll 

360.8 

8.67 

1.158 

2.717 

10091 

3579 

4789 

2.20 

367.2 

6.62 

1.160 

3o270 

10531 

2743 

4620 

2.34 

376.7 

5.27 

1.172 

3.636 

10765 

2360 

4523 

2.42 

381.9 

4.64 

1.177 

4.125 

11028 

1977 

4408 

2.51 

387.7 

3.99 

1.182 

4.635 

11257 

1680 

4300 

2.60 

392.8 

3.48 

1.186 

5.746 

11650 

1249 

4106 

2.75 

401.6 

2.71 

1.193 

7.  160 

12017 

925 

3911 

2.90 

410.0 

2.11 

1.200 

10.482 

12583 

553 

3587 

3.17 

423.0 

1.38 

1.209 

16.422 

13154 

304 

3230 

3.48 

436.3 

0.84 

1.219 

26.911 

13688 

158 

2866 

3.83 

448.9 

0.49 

1.230 

34.625 

13929 

113 

2692 

4.01 

454.6 

0.38 

1.236 

41.552 

14091 

89 

2571 

4.  15 

458.4 

0.31 

1.241 

7* 


TABLE  III- 12 

EXACT  FLOW  VALUES  FOR  RUN  2-C 

Area  Velocity  Pressure     Temp  Mach  Specific  Density  Gamma 

Ratio  (ft/sec)  (lb/ft2)     (Deg  R)  Number  Impulse  x  103  (Frozen) 

(lb- sec/lb)   (lb/ft3) 

1.123  6715  15707   5650    1.33  306.0  24.05  1.530 

1.295  7606  11514   5425    1.55  318.6  18.43  1.155 

1.473  8195  9102   5272    1.69  328.6  15.04  1.157 

1.722  8783  7011   5105    1.84  339.5  12.00  1.160 

2.046  9331  5359   4941    1.99  350.4  9.51  1.146 

2,314  9674  4468   4846    2.09  357,6  8.11  1.136 

2.621  9990  3743   4763    2.18  364.5  6.93  1.127 

3.062  10353  3017   4670    2.29  372»7  5.72  1.123 

3.445  10610  2571   4602    2.37  378.6  4.97  1.124 

4.073  10966  2035   4459    2.49  386.5  4.06  1.185 

4.764  11278  1634   4311    2.60  393.4  3.38  1.189 

6.093  11725  1163   4086    2.77  403.5  2.54  1.196 

9=425  12414  643   3703    3.08  419.2  1.55  1.207 

13.424  12894  400   3412    3.33  430, 4  1.05  1.214 

18,236  13263  266   3173    3.54  439 . 0  0.75  1.221 

29.216  13761  143   2830    3.88  450.7  0.45  1.231 

37.160  13988  104   2664    4.05  456.1  0.35  1.237 

43.092  14119  85   2565    4.16  459.2  0.30  1.241 
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TABLE  III- 13 
EXACT  FLOW  VALUES  FOR  RUN  3-S 

Area         Velocity  Pressure      Temp       Mach         Specific  Density   Gamma 
Ratio         (ft/sec)    (lb/ft2)     (Deg  R)    Number        Impulse  x  103      (Frozen) 
(lb- sec/lb)   (lb/ft3)    , 


1„083  6404  17277  5699  1.27  302.1  26.17  1.342 

1.212  7262  13027  .5466  1.47  312.8  20.61  1.169 

1.444  8127  9319  5230  1.70  326.7  15.46  1.167 

2.008  9289  5415  4849  2.00  348.6  9.73  1.179 

2,625  10020  3605  4562  2.22  363.6  6.90  1.188 

30  194  10481  2704  4361  2.37  373.3  5.42  1.194 

4o  163  11028  1853  4101  2.57  385.3  3.95  1.201 

6o714  11855  954  3665  2.92  403.8  2.28  1.211 

13.303  12787  379  3109  3.41  425.2  1.07  1.225 

22.980  13382  183  2712  3.80  439.0  0.59  1.238 

35.383  13777  103  2425  4.12  448.3  0.37  1.250 
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TABLE  HI- 14 


EXACT  FLOW  VALUES  FOR  RUN  3-CM 

Area         Velocity  Pressure     Temp       Mach          Specific  Density  Gamma 
Ratio         (ft/sec)    lb/ft2)      (Deg  R)     Number        Impulse  x  103  (Frozen) 
„ . _ (lb- sec/lb)  (lb/ft3)  _____ 


1.090 


1.281 


1.401 


1.780 


2.024 


2.412 


6404 


17277        5699 


7575         11605        5365 


7998 


8883 


9273 


9748 


2.781        10100 


3.275        10465 


3.676        10720 


4.023        10910 


4.612        11192 


5.406        11501 


8.911        12335 


15.692        13096 


22.139        13488 


27.270        13704 


32.679        13879 


38.576        14032 


41.708        14100 


9821        5253 


6631        5042 


5472        4958 


4254        4857 


3486        4776 


2801        4711 


2386        4623 


2108       4558 


1737 

1389 

698 

326 

206 

157 

123 

99 

89 


4419 
4255 
3/779 
3301 
3036 
2884 
2756 
2643 
2590 


1.27 


1.54 


1.65 


1.88 


1.98 


2.11 


2.21 


2.31 


2.39 


2.45 


2.55 


2.67 


3.03 


3.43 


3.68 


3.83 


3.96 


4.08 


4.14 


302.1 


317,3 


324.4 


341.1 


349.2 


359.5 


367.3 


375.8 


381.5 


385.8 


392.0 


398.8 


417.6 


435.2 


444.4 


449.4 


453.6 


457.2 


458.8 


26.17  1.286 

18.71  1.168 

16.20  1.154 

11.48  1.126 

9.67  1.120 

7.72  1.112 

6.46  1.122 

5.30  1.161 

4.61  1.143 

4.14  1.048 

3.52  1.203 

2.92  1.205 

1.65  1.212 

0.88  1.221 

0.61  1.227 

0.49  1.232 

0.40  1.236 

0.34  1.240 

0.31  1.242 
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